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Abstract 

Double beta processes in 64 Zn, 70 Zn, 180 W, and 186 W have been searched for with the 
help of large volume (0.1 — 0.7 kg) low background ZnWO"4 crystal scintillators at the 
Gran Sasso National Laboratories of the INFN. Total time of measurements exceeds 10 
thousands hours. New improved half-life limits on double electron capture and electron 
capture with positron emission in 64 Zn have been set, in particular (all the limits are at 
90% C.L.): If]** > 1.1 x 10 20 yr, T*jf + > 7.0 x 10 20 yr, and T°jf + > 4.3 x 10 20 yr. 

The different modes of 2/3 processes in 70 Zn, 180 W, and 186 W have been restricted at the 
level of 10 17 - 10 20 yr. 



PACS: 29.40.Mc, 23.40.-s 
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1 Introduction 

Neutrinoless (Ov) double beta (2/3) decay is one of the low-energy effects which are forbidden in 
the Standard Model (SM) because of violation of the lepton number on two units. However, it is 
naturally expected in many SM extensions. 0z/2/3 experiments offer complementary information 
to those given by neutrino oscillation experiments. While oscillation experiments are sensitive 
to the neutrino mass differences, the 0v2j3 decay rate could determine an absolute scale of 
neutrino mass, and the neutrino mass hierarchy. Besides, double beta decay experiments allow 
to prove the nature of neutrino (that is, if it is a Majorana particle, v = V, or a Dirac particle, 

1 Corresponding author. E-mail address: rita.bernabei@roma2.infn.it (R. Bernabei) 
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v ^ v), existence of right-handed admixtures in weak interaction, and to test a range of other 
effects beyond the SM [H El El SI El El El E|- 

Experimental investigations in this field are concentrated mostly on 2(3~ decays, processes 
with emission of two electrons. Developments in the experimental techniques during the last 
two decades lead to observation of the two neutrino {2u) 2(3~ decay in 10 isotopes with half-lives 
in the range of 10 18 — 10 21 yr, and to impressive improvement of sensitivity to the neutrinoless 
mode of 2(3~ decay up to 10 23 - 10 25 yr [1] . 

Results for double positron decay (2/3 + ), electron capture with positron emission and 
capture of two electrons from atomic shells (2e) are much more modest. The most sensitive 
experiments give the limits on the 2e, e(3 + and 2(3 + processes on the level of 10 17 — 10 21 yr 
[TJ. Reasons for such a situation are: (1) lower energy releases in 2e, e(3 + and 2(3 + processes 
in comparison with those in 2f3~ decay, that results in higher expected T1/2 values, as well 
as provides difficulties to suppress background; (2) usually lower natural abundances of 2f3 + 
isotopes (which are typically lower than 1% with only few exceptions). Nevertheless, studies of 
e(3 + and 2e decays could help to distinguish the mechanism of neutrinoless 2(3 decay (is it due 
to non-zero neutrino mass or to the right-handed admixtures in weak interactions) [10J. 

64 Zn is one of a few exceptions among 2f3 + nuclei having big natural isotopic abundance (see 
Table 1 where properties of potentially 2(3 active nuclides present in zinc tungstate (ZnW0 4 ) 
crystals are listed). This feature allows to build a large scale experiment without expensive 
isotopical enrichment. With the mass difference between 64 Zn and 64 Ni nuclei of 1095.7 keV 
double electron capture and electron capture with emission of positron are energetically 
allowed. 



Table 1: Potentially 2(3 active nuclides present in ZnW0 4 crystals. 



Transition 


Energy 
release 
(keV) (I!] 


Isotopic 
abundance 

(%) m 


Decay 
channels 


Number of nuclei 
in 100 g of ZnW0 4 
crystal 


64 Zn^ 64 Ni 
70 Zn-> 70 Ge 

180 W ^180 Hf 
186 W ^186 0s 


1095.7(0.7) 
998.5(2.2) 
144(4) 
489.9(1.4) 


48.268(0.321) 
0.631(0.009) 
0.12(0.01) 
28.43(0.19) 


2e, e(3 + 
2(3- 
2e 
2(3- 


9.28 x 10 22 
1.21 x 10 21 
2.31 x 10 20 
5.47 x 10 22 



A possible experimental indication of the e(3 + decay of 64 Zn with T^X = (1.1 ± 0.9) x 

10 19 yr was suggested in Ref. [13J. A 07.6 x 7.6 cm Nal(Tl) scintillator and a 25% efficiency 
HP Ge detector, operating in coincidence, were used in that experiment. The excess of ~85 
events in the 511 keV peak was observed with Zinc sample (mass of 350 g, 392 h of exposure 
on the sea level), while no effect was detected without sample or with Copper or Iron blanks. 

During last decade sensitivities of other experiments were not enough to confirm or disprove 
the result of Ref. [13J. CdZnTe semiconductor detector was used to search for 2(3 decays of 
64 Zn in underground measurements performed in the Gran Sasso National Laboratories (LNGS, 
depth of ~3600 m w.e.) during 1117 h in the COBRA experiment; however, small mass of 
the detector (near 3 g) allowed to reach half-life limits at the level of 10 16 — 10 17 yr [14J. 
Low-background experiment with ZnW0 4 crystal scintillator was performed in the Solotvina 
Underground Laboratory (~1000 m w.e.) [15]. Despite the low mass of the detector (4.5 g) 
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and short time of the measurements (429 h), limits at level of 10 18 yr were set for 2(5 processes 
in 64 Zn. 

In the experiment [16J two detectors, HP Ge 456 cm 3 and CsI(Tl) ~400 cm 3 , were used 
in coincidence in measurements with 460 g Zn sample in the underground Cheong Pyung 
Laboratory (^1000 m w.e.). Measurements during 375 h gave the limit on e(3 + decay of 64 Zn 
as: T^ +2y } £l3 > 1.3 x 10 20 yr. Further improvement of sensitivity was reached at the first stage 
[T7] of the experiment presented here. The measurements were performed in the LNGS with 
the help of large ZnWC>4 scintillator (mass of 117 g) over 1902 h. New improved half-life limits 
on 2v and Ou e(3 + decay of 64 Zn were established as: T*jf + > 2.1 x 10 20 yr, T°jf + > 2.2 x 10 20 
yr [IT]. New limits on different modes of double electron capture in 64 Zn have been obtained 
too: Ty% K > 6.2 x 10 18 yr, T^ K > 4.0 x 10 18 yr, and Tffi > 3.4 x 10 18 yr, all at 90% C.L. 

In addition to 64 Zn, ZnW04 scintillators contain a few other potentially 2(3 active isotopes: 
70 Zn, 180 W and 186 W (see Table 1). The best up to date half-life limits on different modes 
and channels of 2/3 processes in these isotopes were obtained in low-background measurements 
carried out in the Solotvina Underground Laboratory (~1000 m w.e.) with ZnW04 and CdW04 
crystal scintillators. The half-life limits of 2(3~ decay of 70 Zn: Ty^ 13 > 1.3 x 10 16 yr, and 
7ffi p ~ > 0.7 x 10 18 yr at 90% C.L. were set in the work [15] . The data accumulated over 
13316 h with 330 g CdW0 4 detector enriched in 116 Cd were analyzed to search for 2f3 decay 
of 186 W at the level of Ty^~ ( 186 W) > 3.7 x 10 18 yr, T°J^~ ( 186 W) > 1.1 x 10 21 yr (all at 
90% C.L.) [18J. Data of 692 h of the measurements with low energy threshold were taken to 
set the following limits on the double electron capture in 180 W: T 1 2 / 1/ 2 2A "( 180 W) > 0.7 x 10 17 yr, 
T°/f £ ( 180 W) > o.9 x 10 17 yr. 

These experiments have demonstrated a good potential of ZnW0 4 (and CdW0 4 ) scintilla- 
tors to search for double beta processes in Zinc and Tungsten isotopes. The main properties of 
ZnW0 4 scintillators are: (i) density equal to 7.8 g/cm 3 ; (ii) light yield ~13% of that of Nal(Tl); 
(iii) refractive index equal to 2.1-2.2; (iv) emission maximum at 480 nm; (v) effective average 
decay time is 24 lis (at room temperature). The material is non- hygroscopic and chemically in- 
ert; the melting point is at 1200 °C. Radiopurity of zinc tungstate crystals has been preliminary 
investigated in [15J. Development of large volume high scintillation properties low-background 
ZnW0 4 crystal scintillators was realized recently pT9l [20] . 

Aim of the present work is to search for 2e capture and e(3 + decay of 64 Zn with the help of 
large low-background ZnW0 4 crystal scintillators. In addition, search for double beta decays 
of 70 Zn, 180 W, and 186 W have been realized as by-products of the experiment. 

2 Measurements 
2.1 Detector 

Two clear, slightly pink colored ZnW04 crystals produced from single crystals grown by the 
Czochralski method were used in our experiment. Two batches of zinc oxide from different 
producers were used to prepare ZnW0 4 compounds for the crystal growth. The crystals used 
for measurements are listed in Table 2. The third scintillator (ZWO-2a) was cut from the large 
0.7 kg crystal (ZWO-2). 
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Table 2: Z11WO4 crystal scintillators used in the present experiments. 



Crystal 


Size (mm) 


Mass (g) 


scintillator 






ZWO-1 


20 x 19 x 40 


117 


ZWO-2 


044 x 55 


699 


ZWO-2a 


044 x 14 


168 



The ZnW0 4 crystals were fixed inside a cavity of 047 x 59 mm in the central part of a 
polystyrene light-guide 66 mm in diameter and 312 mm in length. The cavity was filled up 
with high-pure silicon oil. The light-guide was optically connected on opposite sides by optical 
couplant to two low radioactive EMI9265-B53/FL 3" photomultipliers (PMT). The light-guide 
was wrapped by PTFE reflection tape. 

The detector has been installed deep underground (~3600 m w.e.) in the low background 
DAMA/R&D set-up at the LNGS. It was surrounded by Cu bricks and sealed in a low ra- 
dioactive air-tight Cu box continuously flushed with high purity nitrogen gas (stored deeply 
underground for a long time) to avoid presence of residual environmental Radon. The Cu box 
has been surrounded by a passive shield made of 10 cm of high purity Cu, 15 cm of low ra- 
dioactive lead, 1.5 mm of cadmium and 4/10 cm polyethylene/paraffin to reduce the external 
background. The whole shield has been closed inside a Plexiglas box, also continuously flushed 
by high purity nitrogen gas. 

An event-by-event data acquisition system records the amplitude and the arrival time of 
the events. Moreover, the sum of the signals from the PMTs was also recorded by a 1 GS / s 
8 bit DC270 Transient Digitizer by Acqiris over a time window of 100 /is. To allow a good 
compromise to handle the data files and taking into account the slow scintillation decay of 
ZnW0 4 , 20 MS/s sampling frequency was used during the data taking. 

2.2 Low-background measurements 

The experiment was realized in the LNGS at a depth of 3600 m w.e. The measurements were 
carried out in four runs (see Table 3 for details). The data of Run 1 have been already analyzed 
and published in [17] . 

The energy scale and resolution of the ZnW0 4 detectors have been measured with 22 Na, 
133 Ba, 137 Cs, 228 Th and 241 Am 7 sources. Dependence of energy resolution of the ZnW0 4 
detectors on energy can be fitted by the function: FWHM 7 (keV) = yja + b ■ E y , where E 1 
is the energy of 7 quanta in keV. For instance, the calibration energy spectra accumulated 
with 133 Ba, 137 Cs and 241 Am in the beginning of Run 4 are shown in Fig. 1. The values of 
parameters a and b for the Run 4 are a = 190(40) keV 2 and b = 7.34(35) keV. The energy 
resolutions measured with the ZnW0 4 detectors for 662 keV 7 quanta of 137 Cs are presented 
in Table 3. 

The energy spectra measured with the ZnW0 4 detectors (ZWO-1 and ZWO-2) in the low- 
background set-up are presented in Fig. 2. The spectra are normalized on the mass of the 
crystals and time of the measurements. A few peaks in the spectra can be ascribed to 7 quanta 
of naturally occurred radionuclides 40 K, 214 Bi ( 238 U chain) and 208 T1 ( 232 Th) from materials 
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Table 3: Description of low-background measurements with ZnW0 4 crystal scintillators. Time 
of measurements (t), energy interval of data taking {AE), energy resolutions at 662 keV 7 line 
of 137 Cs (FWHM), and background counting rates in different energy intervals are specified. 



Run 


Crystal 


t 


AE 


FWHM 


Rate (counts/day/keV/kg) 




scintillator 


00 


(MeV) 


(%) 




in energy interval (MeV) 












0.2-0.4 


0.8-1.0 


2.0-2.9 


3.0-4.0 


1 


ZWO-1 


1902 


0.01-1 


11.5 


1.93(3) 


0.27(1) 






2 


ZWO-1 


2906 


0.05-4 


12.6 


1.71(2) 


0.25(1) 


0.0072(7) 


0.0003(1) 


3 


ZWO-2 


2130 


0.05-4 


14.6 


1.07(1) 


0.149(3) 


0.0072(4) 


0.00031(7) 


4 


ZW0-2a 


3292 


0.01-1 


11.0 


1.52(2) 


0.211(7) 







of the set-up (we suppose the PMTs as the main sources of the 7 background). Some part 
of the background counting rate is due to internal radioactive contamination of the crystals. 
Background counting rates in the energy intervals 0.2 — 0.4, 0.8 — 1.0, 2.0 — 2.9, and 3.0 — 4.0 
MeV are presented in Table 3. The background of the detectors will be analyzed in section 3. 

2.3 Mass spectrometric measurements of the impurities in the 
ZnWC>4 crystals 

In order to estimate the presence of naturally occurring radioactive isotopes, the ZnWC>4 crys- 
tals were measured with the help of Inductively Coupled Plasma - Mass Spectrometry (ICP-MS, 
Agilent Technologies model 7500a) at the LNGS. 

Samples of the ZnW04 crystals (ZWO-1 and ZWO-2) were reduced to powder by mechan- 
ical treatment inside a clean polyethylene bag to avoid possible external contamination. The 
samples were etched in a microwave assisted acid digestion technique (Method EPA 3052) using 
nitric acid and nitric-hydrofluoric acid mixtures. The results for both procedures were just a 
partial sample dissolution. The solutions obtained after centrifugation have been analyzed by 
ICP-MS. 

The measurements have been carried out in semiquantitative mode. The instrumentation 
was calibrated with the help of a single solution containing 10 ppb of Li, Co, Y, Ce, and Tl. 
The results of the measurements of the ZnW0 4 samples are presented in Table 4. The accuracy 
of the measurements are on the level of 20-30%. However, the errors for some elements could 
be higher. In particular, in a case of Th interference with tungsten oxide 184 W 16 03 can lead to 
considerable ambiguity. Indeed the measured contamination of the crystal ZWO-1 by thorium 
on the level of 25 ppb corresponds to activity of 232 Th m 0.1 Bq/kg. As it will be shown in 
section 3.1, the 232 Th contamination in the crystal ZWO-1, determined by analysis of data of 
the low-background measurements, is much lower. 

3 Background of Z11WO4 detectors 

Knowledge of the radioactive contamination of the ZnW04 crystals is necessary to describe the 
background in the energy intervals of the 2(3 processes in Zn and W. The time-amplitude anal- 
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Figure 1: Energy spectra measured by ZnW0 4 detector (ZWO-2a) with 137 Cs (main part), 
133 Ba and 241 Am (inserts) 7 sources. 



ysis, the pulse-shape discrimination, and the Monte Carlo simulation were applied in addition 
to the ICP-MS measurements to reconstruct the measured background spectra of the ZnW0 4 
detectors, and to estimate their radioactive contamination. 

3.1 Time- amplitude analysis 

The activity of 228 Th ( 232 Th family) in the ZnW0 4 crystals was determined with the help of the 
time-amplitude analysis^! using the data of Run 2 and Run 3. The arrival time and the energy 
of events were used to select the following fast decay chain in the 232 Th family: 224 Ra (Q a = 
5.79 MeV) -> 220 Rn (Q a = 6.41 MeV, T 1/2 = 55.6 s) -> 216 Po (Q a = 6.91 MeV, T 1/2 = 0.145 
s) — > 212 Pb. a particles in ZnW0 4 scintillators produce a substantially lower scintillation in 
comparison with 7 quanta ((3 particles). Numerically this quenching can be expressed through 
the a/ (3 ratk£l We have used the following energy dependence of the a/ (3 ratio derived from the 
measurements [T5] : a/ [3 = 0.074(16) + 0.0164(40) x E a , where E a is the energy of a particles 
in MeV (E a > 2 MeV). 

Four and seven events of the fast chain 224 Ra— > 220 Rn— > 216 Po — > 212 Pb were found in the data 
of Run 2 and Run 3, respectively. Taking into account the efficiency of the events selection 
(65%) one can calculate the activities of 228 Th in the ZnW0 4 crystals as 5(3) /zBq/kg (ZWO-1) 
and 1.5(6) /iBq/kg (ZWO-2). 

Similarly the limits on activity of 227 Ac ( 235 U family) in the crystals was set by selection of 
the fast decays 219 Rn (Q a = 6.95 MeV, T 1/2 = 3.96 s) -> 215 Po (Q a = 7.53 MeV, T 1/2 = 1.78 

2 The technique of the time-amplitude analysis is described in detail in [21] I22j. 

3 The a/ (3 ratio is defined as ratio of a peak position in the energy scale measured with 7 sources to the 
energy of a particles (E a ). Because 7 quanta interact with detector by (3 particles, we use more convenient 
term "a//3" ratio. 
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Figure 2: (Color online) Energy spectra of ZnW0 4 scintillators measured in the low background 
set-up during Run 2 with ZWO-1 (solid line) and Run 3 with ZWO-2 (dashed line). Width of 
the energy bin is 30 keV. Energies of 7 lines are in keV. 



ms) -> 211 Pb. 

The estimations of radioactive contamination of the ZnW0 4 crystals by 228 Th and 227 Ac 
obtained with the help of the time- amplitude analysis are presented in Table 5. 

3.2 Pulse-shape discrimination between (3{y) and a particles 

As it was demonstrated in [15], the difference in pulse shapes in ZnW0 4 scintillator allows 
to discriminate j(/3) events from those induced by a particles. The optimal filter method 
proposed by E. Gatti and F. De Martini in 1962 [23], developed in [231 125] for CdW0 4 crystal 
scintillators, and successfully used for different scintillation detectors: CeF 3 [26], CaW0 4 [27], 
YAG:Nd [28], ZnW0 4 PSJ, CaF 2 (Eu) [29], PbW0 4 [30], CaMo0 4 [31] was applied for this 
purpose. For each signal /(£), the numerical characteristic of its shape (shape indicator, SI) 
was defined as SI — J2 f(tk) x P{tk)/ J2 f(tk), where the sum is over time channels k, starting 
from the origin of signal and up to 50 /xs, f(tk) is the digitized amplitude (at the time tk) of a 
given signal. The weight function P(t) was defined as: P(t) = {f a (t) — fy(t)}/{f a ('t) + f-y(t)}, 
where f a (t) and / 7 (t) are the reference pulse shapes for a particles and 7 quanta, respectively, 
obtained by summing up shapes of few thousands of 7 or a events. The scatter plot of the 
shape indicator versus energy for the low-background measurements in Run 3 is depicted in 
Fig. 3. The distribution of shape indicators for the events of Run 3 with the energies in the 
energy interval 0.5 — 1.0 MeV is depicted in Inset of Fig. 3. The population of a events is 
clearly separated from events. 

The energy spectra of y(/3) and a events selected with the help of the pulse-shape discrim- 
ination from data of Run 3 are shown in Fig. 4. As it was demonstrated in [15], the energy 
resolution for a particles is considerably worse than that for 7 quanta due to dependence of 
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Table 4: Contamination of ZnW0 4 crystal scintillators measured by ICP-MS analysis. 



Element 


Measured 


Concentration 


Possible 




atomic mass 


of element (ppb) 


interference 






ZWO-1 


ZWO-2 




K 


39 


< 15000 


< 30000 




Ni 


60 


140 


250 




Rb 


85 


1.5 


3.2 




Cd 


111 


670 


11000 




In 


115 


< 2.5 


< 5 




Sm 


147 


< 4 


< 8 




Pt 


195 


16 


10 




Pb 


208 


< io 


< 20 




Th 


232 


25 


< 2 


184 W 16 3 


U 


238 


< 0.5 


< 1 





the a/ 13 ratio on the direction of a particles relative to the ZnW04 crystal axes. Therefore 
we cannot determine activity of U/Th a active daughters in the crystal (equilibrium of U and 
Th chains is supposed to be broken). Fit of the a spectra selected from the data of Run 2 
and Run 3 allows to estimate only limits on activities of U/Th daughters (see Table 5). The 
total internal U/Th a activity in the ZnW0 4 crystals is 0.38(5) mBq/kg (ZWO-1) and 0.18(3) 
mBq/kg (ZWO-2). 

Search for fast decays 214 Bi {Q p = 3.27 MeV, T 1/2 = 19.9 m) -> 214 Po (Q a = 7.83 MeV, 
T1/2 = 164 /is) — > 210 Pb (in equilibrium with 226 Ra from the 238 U chain) was performed with 
the help of pulse-shape analysis of double pulsefl It allows to estimate the activity of 226 Ra in 
the ZnW0 4 crystals at the level of a few /iBq/kg (Table 5). 

3.3 Simulation of j3{^f) background 

The most dangerous radionuclides, which could produce background in a ZnW0 4 detector (in- 
cluding daughters of these radionuclides), are listed in Table 5. As it was mentioned in section 
2.2, radioactive contamination of the PMTs contributes to the background too. The distribu- 
tions of the possible background components were simulated with the help of the GEANT4 
package [22]. The initial kinematics of the particles emitted in the decay of nuclei was given 
by an event generator DECAY0 [33] • To estimate the contribution of the /3(^y) active isotopes 
and of the 7 rays from the PMTs, the measured background spectra of the ZnW0 4 detectors 
(Run 2 and Run 3), after rejection of a events by the pulse-shape discrimination, were fitted by 
the model built from the simulated distributions. Activities of U/Th daughters in the crystals 
were restricted taking into account the results of the time-amplitude and pulse-shape analyzes. 
The activities of 40 K, 232 Th and 238 U inside the PMTs were taken from [33]. The peak in the 
spectrum of Run 3 at the energy « 1.1 MeV cannot be explained by contribution from external 
7 rays (the 1120 keV 7 line of 214 Bi is not enough intensive to provide the whole peak area). 
We suppose presence in the crystal of 65 Zn (Ti/ 2 = 244.26 d, Qp = 1351.9 keV |35j) which can 

4 The technique of the analysis is described in [29] , 
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Table 5: Radioactive contamination of ZnW0 4 scintillators determined by different methods. 



Chain 


Nuclide 


Activity 


mBq/kg) 






ZWO-1 


ZWO-2 


232 Th 


232 Th 


< 0.11 a 


< 0.1 a 




228 Ra 


< 0.2 b 


< 0.05 6 




228 Xh 


0.005(3) c 


0.0015(6) c 


235 JJ 


227 Ac 


< 0.007 c 


< 0.003 c 


238 JJ 


238 JJ 


< 0.1 a 


< 0.08 a 




230 Th 


< 0.13 a 


< 0.07 a 




226 Ra 


< 0.006 a 


0.002(1) a 




210p o 


< 0.2 a 


< 0.06 a 


Total a activity 




0.38(5) a 


0.18(3) a 




4U K 


<l h 


< 0.4 b 




60 Co 


< 0.05 b 


< 0.1 fc 




65 Zn 


< 0.8 b 


0.5 b 




87 Rb 


1.5 d , < 2.6 b 


3.2 d , < 2.3 b 




90g r _90 Y 


< 0.6 b 


< 0.4 6 




137 Cs 


< 0.3 b 


< 0.05 b 




147 Sm 


< 0.01 a 


< 0.01 a 



a Pulse-shape discrimination (see section 3.2) 
b Fit of background spectra (see section 3.3) 
c Time-amplitude analysis (see section 3.1) 
d 1CP-MS analysis (see section 2.3) 
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Figure 3: Scatter plot of the shape indicator (see text) versus energy for 2130 h background 
exposition with the ZWO-2 crystal scintillator (Run 3). (Inset) Distribution of shape indicators 
for the data selected in the energy interval 0.5 — 1.0 MeV. Fit of the distribution by two Gaussian 
functions is shown by solid lines. 



be produced from 64 Zn by thermal neutrons (the cross section of 64 Zn to thermal neutrons is 
0.76 barn [35]) or/and by cosmogenic activation. There are no other clear peculiarities in the 
spectra which could be ascribed to the internal trace contamination by radioactive nuclides. 
Therefore we can obtain only limits on the activities of 40 K, 60 Co, 87 Rb, 90 Sr- 90 Y, 137 Cs, and 
U/Th daughters (see Table 5). The activity of 65 Zn in the crystal ZWO-2 is at the level of 0.5 
mBq/kg. The results of the fit of the spectra (Run 2 and Run 3) in the energy interval 0.1 — 2.9 
MeV, and the main components of the background are shown in Fig. 5 and Fig. 6, respectively. 

The summary of the measured radioactive contamination of the ZnWC>4 crystal scintillators 
(or limits on their activities) is given in Table 5. 

It should be stressed that external 7 rays contribute the main part of the background of the 
ZnWO"4 detectors. For instance ~ 90% and ~ 70% of the counting rate in the energy interval 
0.5 — 1.3 MeV (important to search for the double beta processes in 64 Zn) are determined 
by external 7 quanta in Run 2 and Run 3, respectively. Therefore one order of magnitude 
suppression of background of a ZnW0 4 detector can be reached by improving of the set-up. 
For this purpose we intend to add two high-purity quartz light-guides 10 cm of length each. 
Further advancement needs R&D of radiopure ZnWC>4 scintillators. It should be stressed, no 
special additional low-radioactive measures were carried out to grow the ZnW0 4 crystals used 
in the present experiment. One could assume that the radiopurity of ZnW0 4 is at some extent 
an own property of this crystal. 
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Figure 4: (Color online) The energy spectrum of (3 particles (7 quanta, solid histogram) and 
a particles (dots) selected by the pulse-shape discrimination from raw data measured with the 
ZWO-2 scintillator during 2130 h (Run 3) in the low background set-up. In the inset, the a 
spectrum is depicted together with the model, which includes a decays from 238 U family. The 
total a activity in the ZnWC>4 crystal is 0.18 mBq/kg. 



4 Results and discussion 

There are no peculiarities in the measured energy spectra of the ZnW04 detectors, which can 
be interpreted as double beta decay of Zinc or Tungsten isotopes. Therefore only lower half-life 
limits can be set according to formula: 

lim Ti/2 = N • 7] • t ■ In 2/ lim 5*, 

where N is the number of potentially 2/3 unstable nuclei, rj is the detection efficiency, t is the 
measuring time, and limS 1 is the number of events of the effect searched for which can be 
excluded at given confidence level (C.L.). 

The response functions of the ZnW0 4 detectors for the 2f3 processes searched for and possible 
background components in each Run with different detectors were simulated with the help of 
the GEANT4 code [32]. The initial kinematics of the particles emitted in the decays was 
generated with the DECAY0 event generator [33J. 

We have used different combinations of data recorded in Runs 1-4 to reach maximal sen- 
sitivity to the double beta processes searched for. For instance, the sum of energy spectra 
measured in Runs 2 and 3 was taken to search for double beta processes in Zn and 0z/2/3 de- 
cay in 70 Zn, while all the accumulated spectra were used to set limits on the 2f3 decay in 180 W, 
186 W, and 2z/2/3 decay in 70 Zn (all these processes have a comparatively low energy release). 
Response functions for the double beta processes searched for, the components of background 
describing internal contamination of the crystals, and external 7 rays from PMTs were built as 
sum of models for each Run. 
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Figure 5: (Color online) Energy spectrum of f3(j) events accumulated in the low-background 
set-up with the ZWO-1 crystal scintillator over 2906 h (Run 2) together with the model of 
the background. The main components of the background are shown: spectrum of 40 K (in- 
ternal contamination), and the contribution from the external 7 quanta from PMTs in these 
experimental conditions. 



4.1 Search for double electron capture and ef3 + decay of Zn 

The expected energy distributions for the 2/3 processes in the ZWO-2 detector (Run 3) in 64 Zn 
are shown in Fig. 7. It should be stressed, the registration efficiencies in the whole distributions 
are at least 99.9% for all the processes. 

Two approaches were used to estimate the value of limS* for the two neutrino mode of 
electron capture with positron emission (2vef3 + ) in 64 Zn. In the first one (the so called la 
approach) statistical uncertainty of the number of events registered in the energy region of 
the expected peak was taken as limS*. 11848 events were observed in the energy interval 
490 — 1150 keV of the spectrum Run 2 + Run 3 (see Fig. 8), which gives limS* = 109 counts. 
Considering the related efficiency in the energy interval rj = 0.82, it gives the half-life limit 
Tyf + > 1.0 x 10 21 yr. 

In the second approach the energy spectrum was fitted in the energy range 380 — 1430 keV. 
Background model was composed of 65 Zn, 90 Sy— 90 Y, 40 K, 60 Co, 137 Cs, U/Th inside crystal, and 
40 K, 232 Th, 238 U in the PMTs. The starting and final energies of fit were varied as 380 — 460 and 
1260 — 1430 keV with the step of 10 keV. The fit in the energy region 430 — 1360 keV was chosen 
as giving the minimal value of x 2 /n.d.f. = 121/80 = 1.51. The fit gives total area of the effect 
—232 ±235 counts which corresponds (in accordance with the Feldman- Cousins procedure [3"6] ) 
to limS = 192(65) counts at 90%(68%) C.L. Taking into account ~100% registration efficiency 
for the whole curve in this case, one can calculate the half-life limit: 

T 2 ; 2 £/3+ ( 64 Zn) > 0.70(2.1) x 10 21 yr at 90%(68%) C.L. 

In the same way the half-life bound on the neutrinoless mode was set as: 
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Figure 6: (Color online) Energy spectrum of f3(j) events accumulated in the low-background 
set-up with the ZWO-2 crystal scintillator over 2130 h (Run 3) together with the model of 
the background. The main components of the background are shown: spectra of 40 K, 90 Sr~ 
90 Y, 65 Zn, and the contribution from the external 7 quanta from PMTs in these experimental 
conditions. 



T?jf ( 64 Zn) > 4.3(5.7) x 10 20 yr at 90%(68%) C.L. 

The energy distributions expected for the 2vef3 + and 0z/£/3 + decay of 64 Zn, excluded at 90% 
C.L., are shown in Fig. 8. 

In case of 0z/2e decay, different particles are emitted: X rays and Auger electrons from 
de-excitations in atomic shells, and 7 quanta and/or conversion electrons from de-excitation 
of daughter nucleus. We suppose here that in the nuclear de-excitation process only one 7 
quantum is emitted that is the most pessimistic scenario from the point of view of registration 
of such an event in a peak of full absorption at the Q2/3 energy. 2K, KL, 2L (and other) modes 
are not energetically resolved in the high energy region because of the finite energy resolution 
of the detector (see Fig. 7a). Fit of the measured spectrum in the energy interval 430 — 1300 
keV gives the following limit: 

T°;f £ ( 64 Zn) > 1.1(2.8) x 10 20 yr at 90%(68%) C.L. 

All the limits obtained in the present work for the double beta processes in 64 Zn, as well as 
the results of the most sensitive previous experiments, are presented in Table 6. 

4.2 Search for 2(3~ decay of 186 W and 70 Zn 

The expected energy spectra for different channels of 2/3 decay of 186 W are shown in Inset of 
Fig. 9. 
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Figure 7: (Color online) Simulated response functions of the ZnWCU scintillator (ZWO-2) for 
two neutrino and neutrinoless 2e (a) and e(3 + (b) decays in 64 Zn. One million of decays was 
simulated for each mode. 



The energy spectrum Run 1 + Run 2 + Run 3 + Run 4 was used to search for 2(3 decays 
of 186 W to the ground state and to the first excited (2f) level of 186 0s (E exc = 137 keV). "One 
sigma" approach gives the same limit T^ 13 > 1.1 x 10 21 yr for both transitions (because of 
practically the same registration efficiencies: r] = 99.2% and rj = 97.6% in the energy interval 
400 — 600 keV for transitions to the ground and to the excited states, respectively) . Fit in the 
energy interval 380 - 620 keV gives the limit at 90%(68%) C.L.: 

T o^(i86 W) g g and 2 +) > 2.1(4.2) x 10 20 yr. 

Fit of the spectra accumulated in Run 3 in the energy interval 190 — 450 keV allows to set 
the limit on the 0u2/3 decay with Majoron emission: 

Tjf M1 ( 186 W, g-s.) > 5.8(8.6) x 10 19 yr. 

The best sensitivity to the two neutrino mode of 2j3 decay of 186 W was achieved with the 
spectrum obtained by sum of Runs 3 and 4 (see Fig. 9). Fit in the energy interval 70 — 450 keV 
gives the number of events 2037 ± 1125 which corresponds to \\mS = 3882(3162). It allows to 
set the half-life limit: 

r **^i86 Wj g ^ > 2.3(2.8) x 10 19 yr. 

Fit of the spectrum (Run 1 + Run 2 + Run 3 + Run 4) in the energy interval 190 — 450 
keV gives the following restriction for the half-life relatively to the 2v2f3 decay of 186 W to the 
first excited level 2f of 186 0s: 

T 2l/ 2 2/3 ( 186 W, 2+) > 1.8(3.6) x 10 20 yr. 
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Figure 8: (Color online) The measured energy spectrum of ZnW04 scintillation crystals (Run 
2 + Run 3) together with the response functions for e[3 + process in 64 Zn excluded at 90% C.L. 
Most important components of the background are shown. The energies of 7 lines are in keV. 



Similarly, the limits on 2/3~ processes in 70 Zn were obtained. 

All the results are listed in Table 6 where the best previous limits are also given for com- 
parison. The limits obtained for neutrinoless modes of 2/3 decay of 186 W are on the level of the 
existing results, while the restrictions for the two neutrino models of 2(3 decay in 186 W, and the 
results for 70 Zn were improved by one order of magnitude. The bound on the 0z/2/5 decay of 
70 Zn with Majoron emission was set for the first time. 

It is important to note that in accordance with theoretical calculations in the pseudo SU(3) 
framework [37J which take into account deformation of the 186 W nucleus (in contrast with the 
standard QRPA method), the 2v2(3 decay rate of 186 W could be strongly suppressed (up to 
being equal to 0). Hence, the energy region of 0^2/3 signal of 186 W will be free of the background 
created by the 2v2(3 events, which can reach this region due to the poor energy resolution of 
the detector. The suppression of the 2v mode would be especially important in the search for 
the Ou decay with Majoron emission, whose distribution is continuous, because in this case the 
0z/2/5Ml events will not be distinguished from the 2v background even with the help of the 
high energy resolution detector. 

4.3 2e capture in 180 W 

To set the limits on the 0^28 process in 180 W, the sum of the background spectra of ZnWC>4 
detectors accumulated in all four Runs was used (the low energy part of the spectrum is shown in 
Fig. 10). The least squares fit of this spectrum in the 70—270 keV energy interval gives 133±256 
counts for the 0u2e peak searched for (x 2 /n.d.f.— 12.8/7=1.83), providing no evidence for the 
effect. These numbers lead to an upper limit of 553(389) counts at 90%(68%) C.L. Taking into 
account registration efficiency for this process close to 1, one can calculate the half-life limit: 
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Figure 9: (Color online) Energy spectrum of background obtained by sum of Runs 3 and 4 
together with fit and main components of the background (external 7 ray from PMT, internal 
40 K, 87 Rb, 90 Y- 90 Sr, 238 U). The simulated response function for 2u2/3 decay of 186 W with the 
half-life T^ 13 = 2.3 x 10 19 excluded at 68% is also shown. (Inset) Simulated response functions 
for different double beta processes in 186 W. The distribution for Oi/2/5 decay to the ground state 
of 186 0s is practically indistinguishable from the peak of 0u2j3 decay to the 2 + excited state of 
186 0s. 



Tf/f £ ( 180 W) > 0.86(1.2) x 10 18 yr. 

The same method gives the restriction for the 2v2K process in 180 W: 

r}f( 180 W) > 6.6(9.4) x 10 17 yr. 

The obtained limits are one order of magnitude better than the previous limits obtained in 
the work [18] by using enriched in 116 Cd cadmium tungstate crystal scintillators. The improve- 
ment is reached thanks to absence of the (3 active 113 Cd {Qp = 320 keV [TT] . Ty 2 = 8.04 x 10 15 
yr |38j ) in ZnWC>4 scintillators providing the main background in the 116 CdW04 detectors at 
low energy. The response functions for 0u2e and 2v2K decay of 180 W corresponding to the 
limits obtained in the work [18] are presented in Fig. 10. 



All the half-life limits on 2/3 decay processes in Zinc and Tungsten obtained in the present 
experiment are summarized in Table 6 where results of the most sensitive previous experiments 
are given for comparison. 

All the obtained half-life bounds are well below the existing theoretical predictions [39, 40J; 
nevertheless most of the limits are near one orders of magnitude higher than those established 
in previous experiments. It should be stressed that only a few nuclei among potentially 2e, 
e/3 + , 2{3 + active isotopes were investigated at the level of sensitivity ~10 21 yr. 

Further improvement of sensitivity can be reached by increasing of mass of the ZnWC>4 
detector, suppression of external background, development of ZnWC>4 scintillators with lower 
level of radioactive contamination. 
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Figure 10: (Color online) Energy spectrum of background obtained by sum of Runs 1 — 4 
together with fit in the energy interval 70—270 keV. The simulated response functions for double 
electron capture in 180 W are shown; the half-lives Ty£ K = 7 x 10 16 yr and Ty% £ = 9 x 10 16 yr 
correspond to the the best previous limits reported in [T8] . 



5 Conclusions 

Low background experiment to search for 2(3 processes in 64 Zn, 70 Zn, 180 W, 186 W was carried 
out over more than 10 thousands hours in the underground Gran Sasso National Laboratories 
of INFN by using large volume (117 g, 168 g, and 699 g) low-background ZnWC>4 crystal 
scintillators. 

The new improved half-life limits on double electron capture and electron capture with 
positron emission in 64 Zn have been set, in particular (all the limits are at 90% C.L.): T®J% e > 

1.1 x 10 20 yr, Tyf + > 0.70 x 10 21 yr, and T°jf + > 4.3 x 10 20 yr. The positive indication 
on the (2u + 0u)e(3 + decay of 64 Zn with T 1/2 - (1.1 ± 0.9) x 10 19 yr suggested in [13J is fully 
discarded by the present experiment. To date only three nuclei ( 40 Ca, 78 Kr, and 106 Cd) were 
studied at the similar level of sensitivity (Ti/ 2 ~ 10 21 yr). However, it is worth noting that the 
theoretical predictions are still higher. 

The half-life limits on the 2(5 processes in 70 Zn, 180 W, and two neutrino mode of 2(5 decay 
in 186 W established in our work on the level of 10 18 — 10 20 yr are one order of magnitude higher 
than those of previous experiments. 

We have found ZnWO"4 crystal scintillators extremely radiopure detectors with typical con- 
tamination at the level of /iBq/kg ( 228 Th and 226 Ra), < 0.06 mBq/kg ( 210 Po), total activity 
(U/Th) 0.2-0.4 mBq/kg, < 0.4 mBq/kg ( 40 K), < 0.05 mBq/kg ( 137 Cs), < 0.4 mBq/kg ( 90 Sr- 
90 Y), < 0.01 mBq/kg ( 147 Sm), and < 3 mBq/kg ( 87 Rb). 

An experiment involving ^10 tons of nonenriched crystals (9 x 10 27 nuclei of 64 Zn) could 
reach the half- life sensitivity ~4 x 10 28 yr (supposing zero background during ten years of 
measurements). Such a sensitivity could contribute to our understanding of the neutrino mass 
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Table 6: Half-life limits on 2(3 processes in Zn and W isotopes. 



Transition 


Decay 


Level 




limit (yr) 






channel 


of daughter 
nucleus 


Present work 
90%(68%) C.L. 


Previous results 
90%(68%) C.L. 


64 Zn- 


+ 64 Ni 


0v2e 


g.s. 


> 1.1(2.8) x 10 20 


> 0.7(1.0) x 10 18 p5] 








g.s. 


> 4.3(5.7) x 10 20 


> 3.4(5.5) x 10 18 [T7] 

> 2.8 x 10 16 

> 2.4(3.6) x 10 18 PS] 

> 1.3 x 10 20 HJ5] 






2ue(3 + 


g.s. 


> 0.70(2.1) x 10 21 


> 2.2(6.1) x 10 20 HZ] 

= (1.1 ±0.9) x 10 19 [13] 

> 4.3(8.9) x 10 18 [15] 

> 1.3 x 10 20 [U3] 

> 2.1(7.4) x 10 20 HZ] 


7U Zn- 


* 7U Ge 


0u2p~ 
2u2(3~ 
Oi/2/3-Ml 


g.s. 
g.s. 
g.s. 


> 1.8(3.0) x 10 19 

> 2.3(4.0) x 10 17 

> 1.0(1.4) x 10 18 


> 0.7(1.4) x 10 18 PS] 

> 1.3(2.1) x 10 16 |j5] 


1»"W- 


^18U Hf 


0u2e 


g.s. 


> 0.86(1.2) x 10 18 


> 0.9(1.3) x 10 17 PS] 






2v2K 


g.s. 


> 6.6(9.4) x 10 17 


> 0.7(0.8) x 10 17 [18] 




^186 0g 


Qu2(3~ 


g.s. 


> 2.1(4.2) x 10 2U 


> 1.1(2.1) x 10 2i PH] 






0u2(3~ 


2+ (137.2 keV) 


> 2.1(4.2) x 10 20 


> 1.1(2.0) x io 21 pg 






Oi/2/TMl 


g.s. 


> 5.8(8.6) x 10 19 


> 1.2(1.4) x 10 20 p8] 






2u2(3~ 


g.s. 


> 2.3(2.8) x 10 19 


> 3.7(5.3) x 10 18 p8] 






2v2(3~ 


2+ (137.2 keV) 


> 1.8(3.6) x 10 20 


> 1.0(1.3) x 10 19 p8] 



mechanism and right-handed currents in neutrinoless processes pO]. As well, two neutrino 
double electron capture should be surely observed: in accordance with theoretical expectations 
[3"9"| HO], T1/2 for 2v2e process is predicted on the level of 10 25 — 10 26 yr. 
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